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Abstract 

^~^ ' For calculating large deformations in finite elasticity, we have proposed a method of successive linear approximation, by 

f^ ■ considering the relative descriptional formulation. In this article we briefly describe this method and we prove the existence 

and uniqueness of weak solutions for boundary value problems for nearly incompressible Mooney-Rivlin materials, that arise in 

each step of the method. 
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■<^ . 1. Introduction 

C^ , The constitutive equation of a solid body is usually expressed relative to a preferred reference configuration 

which exhibits specific material symmetries such as isotropy. The constitutive functions are in general 
nonlinear and linearizations can be used as valid approximation only for small deformations. Therefore, 
the problem for large deformations leads to boundary value problems involving systems of nonlinear partial 
differential equations. 



> 

^D • In order to circumvent the difficulties due to the nonlinearities, we have proposed a new method for 

solving numerically the boundary value problem for large deformations. It is based on a successive linear 
approximation by considering the relative descriptional formulation. Roughly speaking, the constitutive 



(N 



f^ ' equations are calculated at each state which will be regarded as the reference configuration for the next 

state, and assuming that the deformation to the next state is small, the updated constitutive equations can 
be linearized. 
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As examples for the proposed method, numerical simulations were done (see [1], [2]) for two classical 
problems concerning Mooney-Rivlin materials, for which the exact solutions are known, namely, the pure 
C^ I shear of a square and the bending of a retangular block into a circular section. The comparison of the 

numerical results with the exact solutions of these two examples confirms the efficiency of our method. 

In the present paper we consider the mathematical analysis of the boundary value problem obtained 
by linearizing the constitutive equations of nearly incompressible Mooney-Rivlin materials relative to the 
present configuration and prove the existence and uniqueness of weak solutions. 

We organize this paper as follows. In Section 2 we introduce briefiy the notion of relative description 
and we formally deduce in Section 3 the linearization of the constitutive function of a nearly incompressible 
Mooney-Rivlin material. In Section 4 we consider a boundary value problem involving a system of partial 
differential equations, that are obtained by linearizing the constitutive equations of a nearly incompressible 
Mooney-Rivlin material, and which corresponds to one of the steps of the successive linear approximation 
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method. The main resuh of this paper is contained in section 5, where we prove the existence and uniqueness 
of weak solutions of this boundary value problem, by considering its variational formulation. For simplicity, 
we restrict our analysis to the two-dimensional case, but the arguments presented can be extended to three 
dimensions. 



2. Relative description and successive linear approximation 

In this section we introduce the notion of relative description, we formally obtain the linearization of a 
general constitutive equation of a solid body respect to this configuration and we describe the successive 
linear approximation method. 

Let kq be a reference configuration of a solid body B, Bq = K(){B), and let 

be the paramctrization of its deformation. Let Kt be the deformed configuration at time t (which we shall 
always refer as the present time), Bt = Kt{B), and 

F{X,t)^VxXiX,t) 

be the deformation gradient with respect to the configuration kq. 

Let Kt be the deformed configuration at time t > t. We define the relative deformation from Kt to k^ as 
the function Xt : Bt —^ Br given by 

Xt{x,T):^X{X,T), xeBt (2.1) 

and the corresponding relative displacement as 

Ut{x,T):^Xt{x,T)-x. (2.2) 

Taking the gradient relative to x in both sides of (2.2), we obtain 

Ht{x,T)^Ft{x,T)-I, (2.3) 

where / is the identity tensor and 

Ht{x,T) := VxUt(a;,T), Ft{x,T) -.^ \/ ^Xt{x , t) 

are called the displacement gradient and the deformation gradient in the relative description, relative to the 
present configuration. 

On the other hand, taking the gradient relative to X in both sides of (2.2), we obtain from (2.1) and the 
chain rule, 

Ht{x, t)F{X, t) = F{X, t) - F{X, t), 

from which we get easily 

FiX,T) = {l + Ht{x,T))F{X,t). (2.4) 
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We can represent this situation by the fohowing picture: 




^^^^^ I+H{t) 





where and hereafter, for simplicity, sometimes only the time dependence is indicated. Position dependence 
is usually self-evident and will be indicated for clarity only if necessary. 

By considering the time t = t + At for small enough At, we can assume that the relative displacement 
gradient is small, 

H{t) := Htix,T), ||if(T)||«l 

Let T be the Cauchy stress tensor given by the constitutive equation 

r = J-.„(F). (2.5) 

Assuming that the operator J^^o is differentiable, we can calculate the linearization of the constitutive 
equation (2.5) relative to the current configuration Kt, and assuming that At is small enough, we have 
formally 

r(T) « Tit) + dJ-«„(^(t))[^(r) - Fit)] = Tit) + d^^,iFit))[HiT)Fit)], (2.6) 

where dJ^Ko(i^) denotes the Frcchet-differential of J^kq calculated at F. For convenience, we shall write (2.6) 
as 

r(r) = T(t) + L(F(t))[if(r)], (2.7) 

where LiF)[H] := dTK„iF)[HF] defines a fourth order elasticity tensor LiF) relative to the current config- 
uration Kt- 

The successive linear approximation method is the discrete construction of the parametrization X(X, t) 
based on the previous arguments. More precisely, let to < • • • < t„_i < t„ < t„+i < • • • be a sequence 
of steps with small enough constant spacing At, where at which step we set t — t„ and t — t„+i. Let 
the deformation gradient FiX,tn) and the elastic Cauchy stress tensor TiX,tn), relative to the preferred 
configuration kjq, assumed to be known. If in any way we calculate the relative displacement Ut^ix,tn+i), 
X G Kf^iB), it allows us the update the new reference configuration Kt„_^-^ relative to the next step by using 
(2.1) and (2.2), i.e., 

X(X,t„+i) := ut„ix,tn+i) + x, X e Kt^iB), 

while the deformation gradient (2.4) and the Cauchy stress (2.7), relative to the preferred configuration Kt^, 
can be determined at instant t„_|_i respectively by 

i^(X,t„+i) :- {l + Ht^ix,tn+i))FiX,tn), 
TiX, t„+i) := T(X, t„) + LiFit„))[Ht„ (a;, t„+i)]. 

Therefore, after updating the boundary data and the eventual body forces acting on the body, we repeat 
the cycle from the updated reference configuration Kt^^-^ ■ 

We remark that this method can easily be extended to constitutive equation T = FiF, F) for viscoelastic 
sohd bodies in general [3]. 
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3. Application to nearly incompressible Mooney-Rivlin materials 

From now on we consider a Mooney-Rivlin material whose constitutive equation relative to the preferred 
reference configuration kq is given by 

T = T^„{F) = -pi + F{F), F{F) - s^B + s^B-', 

where B — FF'^ is the left Cauchy-Green strain tensor and the material parameters si and S2 are constant 
satisfying 

si > and S2 < si. (3-1) 

Remark 3.1: It is usually assumed S2 < < si, the so-called E- inequalities (see [4] and [5]), based on the 
assumption that the free energy function is positive definite for any deformation. Liu [6] has pointed out 
that "any" deformation is unrealistic from physical point of view, and a thcrmodynamical stability anaysis 
only requires S2 < si. Therefore, we shall include the case < S2 < si in our analysis. 

A direct calculation of the Frechet-differential of J^ at F gives 

d:F{F)[H] ^ si{HB + BH^) - S2{B-^ H + H^ B-^) . 

For compressible body in general, the pressure p may depend on the deformation gradient F . However, 
for compressible elastic bodies, we shall assume that the pressure depends only on the determinant of the 
deformation gradient or, by the mass balance, depends only on the mass density p, 



detF(i)' 
where po denotes the mass density in the referential configuration kq. 

For time t — t + At and from (2.4), we have 

p(r) - p{t) = po(dctF(r)-i - dct F{t)-') = pit) (dct F(t)F{T)-' - l) 
= p(t)(det(/ + i7(T))-i -1) = -p{t)tYH{T)+o{2), 
where tr H means the trace of H and o(2) denotes higher order terms in the small displacement gradient 
H{t). Therefore, assuming that p is differentiable as function of p, we have 

P(T)-P(t)=(^) (p(r)-p(t))+o(2) = -(^p^) tri/(r)+o(2), 
or 

p{T)=p{t)-P{t)tTH{T)+0{2), 

where /3(i) = p{t){dp/dp)t is a material parameter depending on the mass density p. 

A body is called nearly incompressible if its density is nearly insensitive to change of pressure. Hence, if 
we regard the density as a function of pressure, p = p{p) , then its derivative with respect to the pressure is 
nearly zero. This means that, for nearly incompressible materials, the parameter (3 must be large, 

/3==/3(a;,t)>l, Vx e Bf 

Therefore, the Cauchy stress tensor relative to the current configuration Kt is given by 

r(r)=T(t) + L(F(t))[i7(r)]+o(2), 

where 

L{F)[H] = P{tT H)I + si{HB + BH'^) - S2{B-^H + H^ B'^) 

and the first Piola-Kirchhoff stress tensor at time r relative to the current configuration nt is given by 

T,,{t) ^ dot Ft{T)T{T)Ft{Ty^ ^det{I + H)T{T){I + H)-^ 

= [1 + tri7 + o(2)] [T{t) + L{F{t))[H] + o(2)] [l - H'^ + o(2)] (3.2) 

= T{t) + {tr H)T{t) - T{t)H^ + L{F{t))[H] + o(2). 



4. Linearized boundary value problem and its variational formulation 

For simplicity, we denote by k the current configuration kj, 17 = '^('S) be the bounded domain of M"^ 
representing the interior of the region occupied by the body at current configuration k at the present time 
t, Tq = T(t) and Bq = B{t). Let dVl = Fi U r2 U F3, n^ be the exterior unit normal to dO. and g be the 
gravitational force (per unit mass). 

We consider the following boundary value problem for the relative displacement u = u(x,t), 

- divr^(T) = p{T)g in fi X R, 
T^{T)n^^f{T) onFi, 

u(T)-n, =0 onF2, (4.1) 

T^iT)n^, X n^ = on F2, 
u{t) = on F3, 

where div is the divergence operator with respect to x, Ti^{t) = Tf^{x, r) is the Piolla-KirchhofF stress tensor 
at time r relative to configuration k at the present time t, which, up to linear terms in relative displacement 
gradient H = H{t) = Wxu(t), is given by (see (3.2)) 

T« = To + (trH)(To + 131) - T^H^ + s^{HBq + BqH'^) - S2{B^'H + H^B^^), 

and /(r) is the surface traction (per unit surface area). 

Remark 4.1: At every time step, the idea of formulating the boundary value problem in the form (4.1) is 
similar to the theory of small deformations superposed on finite deformations (see [7], [5]). In this manner, 
either we are interested in the evolution of solutions with gradually changing boundary conditions resulting 
in large deformation, or, we can treat the boundary values of finite elasticity as the final value of a successive 
small incremental boundary values at each time step (see [8]). 

The boundary value problem (4.1) can be formulated as a variational problem. Indeed, let fJ be a smooth 
enough bounded domain in M? and define the space 

V = {u e {H^i^yf ; M • Wk = on F2 and n = on F3}. 

Taking formally the inner product of both sides of the equation in (4.1) by to G V and integrating over O, 
we obtain after integration by parts, 

/ tr{K[H]W^) dx^ I f{T) -wdV- I tr(ToVF^) dx, 

where we are denoting H = WxU, W = VxW and K[H] is given by 

K[H] := [ti H){To + 131) - TqH^ + si{HBo + BoH^) - S2{B^^H + H^B^^). 
Therefore, for u, w Cz V we consider respectively the bilinear and the linear forms: 

£{u,w) := I tT{K[H]W^)dx, 

J\f{w) := / /(t) -wdT- / tr(ToT4^^) dx + p{T)g ■ w dx. 
Jr, Jn Jn 



(4.2) 



We notice that the forms £ and M can be written in terms of coordinates by 

Jndxk\ -"JOxj Ji2 dxkdxj 

f [ dui , , , , duj \ dwi , 

J^iw)^ / f^WidV- I [Ta]ij-—^dx+ / pgiWidx, 

where in the above formulas we have used the standard summation convention for repeated indices. 
Then, the variational problem is to find the solution u G V such that 

C{u,w):^M{w), yweV. (4.3) 

In order to prove that the solutions of (4.3) is a weak solution of (4.1), the following result concerning 
existence of a normal trace is useful. 

Lemma 4.2: Let il C R^ be a bounded Lipschitz domain. If F e L^(fi)^ satisfies div F e L'^i^), tlien 
F ■ n^ can be defined as an element of H^^^^{dfl) and there exists a constant Ci > depending only on fi 
such that 

\\F-n4H-i;2 <Ci(||F||2 + ||divF||2^ 



Proof: Assume that F e C^{fl) n C°{Ti). Then, using integration by parts, por any ip e C'-iO,) n C"{Q) we 

have 

/ F{x)-\/i:{x)dx+ / div F{x)tlj{x) dx = / tl;{x)F{x) ■ n^{x) dT. 
Jn Jn Jdn 

Therefore, denoting the right hand side of the above identity as {F ■ n^ ; 7o(V')), with the brackets meaning 

the duality between i/^^/^(90) and H^^^{dfl) and 70 : H^{il) -^ H^^^{dfl) being the trace operator, we 

have 

|(F • n, ; 7o(V'))| < IIFII2IIV7MI2 + II divF||2||^||2. 

It is well-known that, for a given ip e H^^^{dn) we may choose ip e H^{Q) such that 7o('0) — V ^nd such 
that ll'011-f/i 1^ C'i||(^||/ji/2, where the constant C'l depends only on fl. Hence, 

\{F ■ n, ; ip)\ < Ci{\\F\\2 + \\ divf II2) ||(^||hi/2. (4.4) 

This means that 

||-F-n«||g-i/2<Ci(||F||2 + ||divf||2). 

When F is no longer in C^{n) C'^(f7), using a density argument (see [9]), we can find a sequence {-FnjneN 
inCi(f7)nCO(f7) such that 

Fn^F in L^(^)^, divf„^divf in L^{n). 

Inequality (4.4) shows that {F„ • n^jngN is a Cauchy sequence in iJ^^/^(951), whose limit, which is inde- 
pendent of the particular choice of the sequence {F„}„gN, wiU be denoted by F • n^. This finishes the proof. 
D 
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Lemma 4.3: Let O C M'^ be a domain of class C^. We assume that (3,pa e L°°{n), p e L^(fi) and 
To, Bq e i°°(r2, M3(R)), wiiere M3(R) denotes the set of 3x3 real matrices. If u is a solution of (4.3), then 
u is a weak solution of (4.1). 

Proof: Let w e V be a solution of (4.3). Then, H := V^^u e L'^(Q,, M3(R)), which imphes that T^ = 
To + K[H] e L^{n, M3(R)). Since C^{nf C V, we have 

- / divT^wdx^ / tr(r„T^'^) dx ^ pg ■ w dx, Vw e C^^i^f, 
Jq. Jn Jn 

where W = V^ w and the partial derivatives in div are taken in the sense of distributions in il. Hence, u 
satisfies 

- div T^ = pg 

in the sense of distributions. Moreover, since we are assuming that p E i^(51), it follows from the density of 
C(f (17) in 2.2(0) that divT^ e L'^{nf. From Lemma 4.2, (4.3) reduces to 



T,n^-wdr^ / fwdT, (4.5) 

where the above surface integral on dfl are taken in the sense of the duality between H^^''^{dn)^ and 
H^^^{dfl)^. In particular, for any w gV such that to = on r2, we have 

/ {T^n^ -f)-wdr = 0, 

which gives the Fi-boundary condition in (4.1). So, (4.5) reduces to 

T^n^ ■wdr = 0, Vw eV. (4.6) 



r. 



In order to show that (4.6) gives the r2-boundary condition in (4.1), let ip e _ffQ(J7), (p < 0, he the first 
eigenfunction of —A and define 

Wo{x) ■.= 'Vx'p{x)\'Vxf{x)\^'^, xen. 

Since fl is of class C^, we can extend Wq to the boundary dfl and we have from the maximum principle that 
100(3;) = ni^(x), for almost all x £ dil. Let w € H^{il)^ be an arbitrary function which vanishes on Fa and 
consider w — Wq x w. Then, it is clear that w d V, since to = on F3 and 

w ■ n^ Ir^ = (n^ x w) ■ n^ \r., = -(n^ x n^) ■ w\r^ == 0. 
Therefore, from (4.6), 

= / Ti.,n^ ■ wdT = / T^n^ ■ (n^ x w) dV = / (r^n^ x Uk) ■ w dT 

JT2 '^^2 '^T2 

and the proof is complete. D 



5. Existence and uniqueness of solution in t^vo-dimensions 

Let ri be a bounded Lipschitz domain of K^ whose boundary dil = Fi U r2 U Fa, with meas(Fi) 7^ for 
i = 1,2,3, and consider the space 



V = {u = {ui,U2) e H^{9.Y ; tt • n^ = on F2 and u = on F3}. 
For u,w G V, we introduce 

(u|i;) := / iyui{x) ■Vvi{x) +Vu2{x) ■Vv2{xyf dx, 

\\u\\l:=\\Vui\\l2 + \\Vu2\\ 



(5.1) 



(5.2) 



2 

L2i 



where || \\i^2 is the usual L^-norm. It is well-known that the Poincare inequality holds if nieas(F3) 7^ 0, i.e., 
there exists a constantc C such that 

\\u\\l>C\\u\\l2, VneV. 

In this case, (•; •) and || ||v define an inner product and a norm in V, respectively. 
From now on we assume that 



where by 5*2 (K) we denote the set of all symmetric and positive definite 2x2 matrix, and we set 

To := -pol + siBo + S2B^\ 

It is clear that the forms C and Af defined in (4.2) are continuous in V. 
Recalling that H and W arc 2x2 matrix whose entries are given by 

dui r,,„ dw 

the bilinear form C{u, w) defined in (4.2) can be written as 



(5.3) 



[^^^■'" = 9^' t^l^^'^a^' "'^e^' 



C{u,w)^ I A{x;H{x),W{x))dx, 
Jn 



where 



A{x; H, W) := ti{H) tr (Tq + I3I)W' - ti{ToH^ W 



TTiirT\ 



Si tr 



T\TjrT 



(HBo + BoH' )W 



S2 tr 



Td-1\ti/T 



iBo'H + H'B^')W 



In particular, for W = H we have 

A{x- H, H) = tT{H) tr [(To + l3I)H^] - tr(ToH^H^) 



Si tr 



{B^'H + H^B^')H 



[HBo + BoH^)H^\ - 32 tr 
Hence, to prove that C is coercive, it is sufficient to show that there exists a > such that 

A{x;H,H) >a\\Hf, Va; en, 
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(5.4) 



i.e., it suffices to show that the biiinear form A{x; H, W) is uniformly coercive as function of 2 x 2 matrices. 
Furthermore, a direct calculation (see the Appendix) gives that the coercivity oi A{x,H,W) is equivalent 
to the semipositivity of the matrix A{x) — al, for all x G ft, with A{x) given by 



A{x) = 



/,3+2si7i-2s27r^ ;3+itrTo 

/3+itrTo l3+2sij2-2s2j~ 


\ 







2si tr_Bo-2s2 trB^ -trTo si (72 -7i)-S2(77 -I2 ) 

trTo 



si(72-7l)-S2(7i ^-72 ^) 



(5.5) 



where 71 e 72 are the eigenvalues of Bq and / is the 4x4 identity matrix. Therefore, we can also write 



C{u,u)^ X{x)'^ ■A{x)X{x)dx, 
Jn 



where, following the notation introduced in the Appendix, 



H{x) 



and 



X{xf := {a,c,b,d) = 



dx\ dx2 

du2 du2 

dx\ dx2 ^ 




= 


a b + d 
b-d c 






dui du2 1 


dui du2 
dx2 dxi 


1 


dui 


du2 


dxi ' dx2 ' 


2 


'2 


dx2 


dxi 



We notice that if A{x) — al is uniformly semipositive in 17, then 



X{xf ■ A{x)X{x) > a 



dui 
dxi 


2 
+ 


du2 
dx2 


^ 1 

+ 2 


dui 
dx2 


' 1 

+ 2 


du2 
dxi 



and consequently. 



C{u,u)>^{\\Vu4l2 + \\Vu2\\l2) = ^\\u\\l. 



In order to analyze the matrix (5.5) and in view of the conditions (3.1), we must distinguish two cases: 
S2 < < si and < S2 < si. In both cases, we fix a constant k > max{0, S2S1 } and take e := si — S2k^^- 
Now, let flo = a(){x) and bo = bQ{x) be the functions defined by 



aa 



2s2 tr Bq" ^ - 2 ( s 1 VdctSt^ - S2 V dot 5(7 M , 



60 :=-2.S2trB/ + 2(siVdct~B^-S2VdctBo^ 



(5.6) 



Assuming that deti?o > k, we have the following inequalities 



&o - flo > 4V-S1S2 if S2 < 0, 



&o - flo = 4 si 



S2 



dot Bo 



v/dct Bo > 4eVA: if S2 > 0. 



(5.7) 



Notice that the above inequalities indicate that the interval [ao(a;), foo(a;)] has nonempty interior for all x E il 
if det i?o ^ ^, and this will be essential in proving the next theorem. 

Finally, we set 

d = sup 



/ trSo 



xen \ Vdct Bo 
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Theorem 5.1: Suppose that detBo > k. Let a > such that 



{2eVk ifs2>0, 

and assume that po satisGes the condition 

ao{x) + ad < —2po(x) < bo{x) — ad, V x E il. 



(5.8) 



(5.9) 



Then, there exists a constant j3q = /3o(si,S2,a) > such that the matrix A{x) — al is uniformly positive 
semidehnite in il, provided that /3{x,t) > /3o for almost all a; e 17. 



Remark 5.2: Since 2-\/det Bq < tr Bq, it follows that d > 2. Hence, if a satisfies (5.8), we have necessarily 

V-S1S2 if S2 < 0, .g -^Qx 



a < , ^ 

. eVk if S2 > 0. 

Proof of Theorem 5.1: The nonzero entries of the matrix A are 



All := ,3 + 2si7i - 2s27i \ 



A, 



13 + 2si72 - 2s272 



-1 



A12 :=^+-trro, 



== 2si tr Bq - 2s2 tr Bg"^ - tr Tq, 

= trTo, 

= si(72 - 71) - S2(7r^ - 72^^), 



where 71 e 72 are the eigenvalues of Bq. To simplify the notation, we introduce the functions /, 5 : (0, +00) — > 
M, as 

7(7) := S17- 527^^ and 5(7) := S17 + S27~^ 

A necessary and sufRcient condition for the matrix A — al he positive semidehnite is 

minJAii - a, (An - a){A22 - a) ~ AI2, A33 - a, {A33 - a){Ai4 - a) - Aj^^ > 0, Va; e n. (5.11) 

It is clear that the condition (5.11) implies, in particular, A22 — a > and ^144 — a > 0, since A is symmetric. 
Step 1: Analysis of the first block of A: 



In the case S2 < 0, we have 7(7) > \/—siS2 for all 7 > 0. So, 



Aii-a^ l3-a + 2/(71) > /3 - a + 2a/-siS2 > /S - a. 

In the case S2 > 0, we can assume without loss of generality that 71 > 72. Then, as si — S2/ dot Bq > fc, we 
have 

S2 



/(71) > si7i - — > fc7i > 0, 
72 



and 



All - a == /3 - a + 2/(71) > /3 - a + 2/071 > /3 - a. 
Therefore, in the two cases. An — a > ii /3 > a. 

On the other hand, if we denote fi — /(7i) and gi — g{ji), i — 1,2, we get 

1 



l3+-tTTQ 



[All - a){A22 -a)-Al2^W-a + 2/i)(/3 -a + 2/2) 

= iP~a) [2(/i + /2) - 2ia - pq) - {gi + 52)] + 4/1/2 - 
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{a - Pq) 



.91 + .92 



n 2 



Since 

2(/i + /2) - (.gi + 92) = Si trBo - 3s2 ir B^\ 

,91+92 1 . rp 

we have 

(All - a){A22 -a)- Al^ - (/3 - a) [si tr Bo - 3s2 tr B^^ - 2{a - po)] + 4/1/2 - 
Therefore, if 



1 ^2 

a + - tr To 



-2po < -2Q; + sitrBo -3s2trS(7\ Va; G f7, (5.12) 

it foUows that {An — a){A22 — a) — ^12 > for /3 > large enough. 
Step 2: Analysis of tie second block of A: 
By the definition of To, we have 

A33 - oi = si tr Bq - 3s2 tr Bq^ + 2po - a. 

Hence, A33 — a > if, and only if, 

-2po<-a + sitrBo-'is2tTB^\ (5.13) 

We notice that A44 — a = —2poSi tr Bq + S2 trSo" — a, which implies that A44 — a > if, and only if, 

-2po > a-sitrBo- S2trBo"^ (5.14) 

The conditions (5.13) and (5.14) can be expressed by 

— si tr Bq — S2 tr Bq + a < Q < —a + sitr Bq — 3s2 tr Bq , 

It is noteworthy that (5.12) implies (5.13). Moreover, the interval 

[-si tr Bq - 32 tr Bq^ + a , -2a + si tr Bq - 3s2 tr Bq^] (5.15) 

is not empty, because if we denote 

{a* := -si tr Bq - S2 tr Bq^ , 
6* := si tr Bq — 3s2 tr Bq^ , 

it follows that 

5* — a* — 3a = 2si tr Bq — 2s2 tr Bq — 3a. 

and hence, from (5.10): 



1) in the case S2 < we have 6* — a* — 3a = 2/i + 2/2 — 3a > A^/ —siS2 — 3a > 0, 

2) in the case S2 > we have 6* — a, — 3a = 2 tr Bq I si ) — 3a > AeVk — 3a > 0, 

V det Bq J 

which implies that the interval defined by (5.15) is not empty. 
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Now, according to the notation introduced above, we have 

^33 := /i + /2 - 2s2 tr Bq-i + 2po 
Am -=91+92- 2po 
A34 := .92 - 51 



So, 



{A33-a){A4,i -a)-Al 



fi+f2-2s2trBo^+2po-a 
[F + 2po][G-2po]-[92-9i? 



91 + 92 -2pQ - a 



[92 - .91 ] 



where we are denoting 



Hence, 



F := /i + /2 - 2s2 tr Bq -a and G :== 51 + 52 - «• 



(A33 - a) (^44 - a) - ^34 =FG + 2po{G - F) - (32 - 5i)' - 4pg. 
For X := — 2po, we have (A33 — a)(yl44 — a) — A34 > if, and only if, 

X^-(F- G)X + (32 - .91)^ -FG<0. 

The above incquahty holds if the discriminant of the binomial (5.16) is positive. In fact, we have 

(F - Gf - 4[(.g2 - 9if - FG] = [,h + /2 - 2^2 tr B^' + 51 + 32 - 2a] ' - 4(32 - 5i)' 
= [2si trBo - 2s2 trBp-i - 2a]' - 4[si(72 - 71) + ^2(72"' - 7r')]' 
= 4[si(7i + 72) - S2(7f ^ + 72^^) -a] - 4[si(72 - 71) + 32(72"^ - li^)] 
= 4[2si72 - 2s27r^ ^ Q^] [2si7i - 2s272"^ ^ Q^] 
= 16[si72 - S27j7^ - "/2] [si7i - 5272^^ ~ Q^/2] 
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/(VdetSoj --{sitTBo-S2trB-') + — 



Note that 



si trSo — S2 trSQ =sitri?( 



- S2 



detSo 



tr Bo si 



S2 



dctB, 



(;^)^<^^)^ 



To simplify the notation, consider 



C -.^ f{y/dctBo) and D := 



trBo 
Vdet Bo ' 



Then, 



C^ 



2 4 



-CD == CM 1 



2C 



Note also that, from (5.8) we have: 



l)ifs2<0, aD <ad< 2y/-siS2 < 2f{^jActBo) = 2C; 



2) if S2 > 0, aD <ad< 2eVk < 2VdctSo si - , "' ) = 2f(JdctBo) = 2C, 

' dot Bo ' 



S2 



(5.16) 
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which implies that, in both cases, < aD/2C < 1. So, by calculating the roots a„ and 6„ of the binomial 
(5.16), we get, 



-2s2tr 5(7^ + 2 



-2s2trBr;^ -2 



/(VdctBo)] -^(sitrSo-S2trB, 



} 



T 



(5.17) 



/(VdetSo)] -^(sitrBo-satrB, 



^ 



T 



and the condition (^33 — a)(A44 — a) — A34 > is equivalent to a^ < — 2po < &o 
We can rewrite (5.17) as 



-2s,tTB^' + 2\IC^-^DC+^, 



a„ = -2s2 tr Bg-i - 2\/C^ - -DC + — , 



so that 



C^~^^CD + ^>JC^-^^CD^C^1-^^§>C-^^D. 



Therefore, from (5.8) and (5.6), we have 

ba ^ bg — ad and Aq < Oq + o^d,- 

Notice also that, from (5.7) and (5.8), we have ba — aa > bo — ao — 2ad > 0. Thus, to conclude the proof, it 
suffices to show that, under the hypothesis (5.8), the following inequalities hold: 

Qt: + a < qq + ad and 60 — ad < 6* — 2a. 
Indeed, first note that 



2Vdet~Bo < trBo and 2y det B(7^ < tr 5(7^ 
Therefore, in the case S2 < 0, we have 



sitrSo - S2trB(7^ > 2siv/detBo - 2s2y dot B^^ = 2f{^/dctBa) 
from which we conclude that 



si tr So - 3s2 tr B(, ^ > -2s2 tiB„^ + 2f{^/dctBa) 
and so, &, > fco- On the other hand, it is easy to see that 



-sitrBo + S2trBo"^ < -2si Vdct Bq + 2S2J Act B^^ = -2f{y/dctB^), 



which implies that, 
and so, a, < oq. 



-sitr^o - S2tTB~^ < -2s2trB(7^ -2/(v/detSo 
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In the case S2 > 0, we have 



S2 

tiBo 



6, - 2a > 60 - ad <^=^ ti Bq si - , ^ - 2a > 2/(VdctA^) - arf 

* net Bo ' ^ ' 



/(VdetSo) - 2a > 2/(\/detBo) - ad 



Vdct So 



Since trSg > 2Vdct Sq and d > 2, it follows that 6* — 2a > 60 ^ Q^c' holds for all a > 0. 
Likewise, 



a < ao + ad -^==^ -si tr So + S2 tr i?Q ^ + a < — 2/(-\/dct So) + ad 

^-trSo)-..(^-trV)<(d-l)a 

■^ V ) (2\/detBo - tr Bo) < (^ - l)a 



si 



detS, 
and a* + a < flo + ad holds for all a > 0. This finishes the proof. D 

Remark 5.3: The above considerations permit us to conclude (by Lax-Milgram Lemma) that the boundary 
value problem (4.1) admits a unique solution. In fact, 

Corollary 5.4: Under the hypothesis of Theorem 5.1, the variational problem (4.3) admits a unique solution 
ueV. D 

Remark 5.5: Theorem 5.1 gives a sufficient condition for the existence of a unique weak solution of the 
boundary value problem (4.1) corresponding to each time step of the successive approximation. As we are 
supposing that the material is nearly incompressible, it is reasonable to expect that det Bq ~ 1- This implies 
that, if 71 and 72 are the eigenvalue of Bq, 71 ~ I/72 and tr Bq^ « 71 + I/71. So, the hypothesis (5.9) does 
not means that we are assuming that po is small. On the other hand, numerical experiments show that the 
hypothesis (5.9) can be very restrictive in the presence of gravitational body forces. In this case, we can 
incorporate the potential of the gravitational force into the pressure, and analyze the re-formulated problem. 

Remark 5.6: The previous results hold if we assume that Fa = 0. In fact, unlike the space V introduced in 
(5.1), we must consider 

V= {ue {H^{n)f ; u • n« = on Ta}. (5.18) 

However, in this case, it is necessary to assume that the domain fl satisfies a geometric property to ensure 
that (5.2) is a norm. This can be done by supposing that 17 has the following property: There is no constant 
vector c G R^ such that c ■ ni^(x) = 0, Va; € Fa. 

6. Appendix 

Without loss of generality, we can assume that Bq is a diagonal matrix, given by 

'71 



To = 



\ 72 

and, in this case, 

'h 0\ /-po + /(7i) 

t2y ^ -po + fM 
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where f{j) = Si7 — S27 "'"• 

Writing the quadratic form (5.4) as 

A{x,H, H) ^ Ai{x, H, H) + A2{x, H, H) + Az{x, H, H) + Ai{x, H, H), 

where 

' Ai{x, H, H) = tr(i7) tr [(Tq + PI)H^] , 

A2{x,H,H) = -tT{nH^H^), 

Asix, H, H) ^ si tr UhBo + BnH'^)H^] , 

Ai{x,H,H) = -S2ti\{Bo^H + H^Bo^)H^ . 
and the matrix H a.s H = E + R, where E ^ ^{H + H^) and R = ^{H - H'^), with 

R 



E ■ 



a b 

b c 



d 
-d 



we obtain, 

1) tr{ToH'^ + /3H'^) = tr{ToE + I3E) = ah + cfa + I3(a + c), which gives 

Ai{x,H,H) = {a + c){ati + ct2) + I3{a + cf . 



(6.1) 



2) Since TqH^ H'^ = To(E^ + R^) - Tq{ER + RE), we have tr(Toi/^i/^) = ii[TQ{E'^ + R^)] and a direct 
calculation gives 

A2{x,H, H) ^ -ti{a^ + b^ - d^) - t2ib^ + c^ - d^). 



(6.2) 



3) We notice that 

J BoH^H^ = Bo{E^ + R^)- Bo{ER + RE), 

\ HBqH'^ ^{EBoE-RBoR) + iRBoE-EBoR). 
Since Bq{ER + RE) and RBqE — EBqR are skew symmetric, we have 

trUnBo + BoH^)H^] = tT[BoiE^ + R^) + (EBqE - RBqR)] , 
and a direct calculation gives 

A3{x,H,H) = 2si 7ia2 + 72c2 + (7i+72)&^ + (72-7i)&c? . 

4) As before, 

J B^^HH'^ = B^\E^ - R^) + Bo{RE -ER + RE), 

\ H'^BqH^ = {EB^^E + RB^^R) - {EB^^R + RB^^E), 
a direct calculation gives 

Ai{x,H,H) = -2s2[l^^a'+l2^^ + {l^^ + l2^)b' + {l^^ - l2^)bd]. 



(6.3) 



(6.4) 



Therefore, by denoting X ~ {a,c,b,d)'^ and considering (6.1)-(6.4), we can express the quadratic form 
(5.4) as 

A{x,H,H) = X'^ ■ A{x)X, 

where A(x) is the matrix 



A{x) = 



'/3+2si7i-2s27r^ /3+itrTo \ 

/3+itrTo /3+2si72-2s272"^ 

2sitrBo-2s2trS-i-trTo si (72 -7l)-S2 (7r^ -72^^) 

si(72-7i)-S2(7j7"'-72^"') trTh / 
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